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Abstract 

Polynomial and spline quasi-interpolants (QIs) are practical and effective 
approximation operators. Among their remarkable properties, let us cite for 
example: good shape properties, easy computation and evaluation (no linear 
system to solve), uniform boundedness independently of the degree (poly- 
nomials) or of the partition (splines), good approximation order. We shall 
emphasize new results on various types of univariate and multivariate poly- 
nomial or spline QIs, depending on the nature of coefficient functionals, 
which can be differential, discrete or integral. We shall also present some 
applications of QIs to numerical methods. 



1 Introduction 

A quasi-interpolant of / has the general form 

Qf = Ma 

where {B a , a € A} is a family of polynomials or B-splines forming a partition of 
unity, and {/j, a (f),a 6 A} is a family of linear functionals which are local in the 
sense that they only use values of / in some neighbourhood of = supp(B a ). 
The main interest of QIs is that they provide excellent approximants of functions 
without solving any linear system of equations. In the literature, one can find the 
three following types of QIs: 

(i) Differential QIs (abbr. DQIs) : the linear functionals are linear combinations of 
values of derivatives of / at some point in E Q . 
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(ii) Discrete QIs (abbr. dQIs) : the linear functional are linear combinations of 
values of / at some points in the neighbourhood of S a . 

(iii) Integral QIs (abbr. iQIs) : the linear functionals are linear combinations of 
weighted mean values of / in the neighbourhood of S Q . 

We shall present various types of univariate and multivariate polynomial and spline 
QIs, mainly dQIs and iQIs, which were recently introduced in the literature. For 
polynomial QIs, we only present QIs which are close to the orifginal Bernstein or 
Durrmeyer operators (for other types of QIs, see for example [35] [36]). 
The prototype of polynomial dQIs is the classical Bernstein operator 

n 

B nf = Y J fC-)t ) 

i=0 

where {b\ n \x) — C l n x l {l — x) n ~\0 < i < n} is the Bernstein basis of the space 
P„ of polynomials of degree at most n (the C % n are binomial coefficients). 
The prototype of polynomial iQIs is the Durrmeyer operator [33] 

n 
i=0 

where b\ n) = b\ n) / ^ b\ n) = {n+ l)b<f l) and (f,g) = fg. Both can be extended 
to the multivariate case, either on the hypercube or on the simplex. Another 
extension consists in adding a Jacobi weight in the scalar product. 
The prototypes of spline DQIs are de Boor-Fix QIs [11] and their various univariate 
and multivariate extensions 

Qf = J2 X ^ B r 

je.J 

Here {Bj j € J} is a family of univariate B-splines of degree m on a nonuniform 
sequence of knots Assuming that £j = supp(Bj) = [tj- m ,tj+i], we set 

E m = {—m + 1, . . . , 0} and we define ipj (t) = Y\ reEm (tj+ r —t)& P TO for all j 6 J. 
For any t e Sj, the coefficient functionals are 

-. rn — 1 

W) = T^Tv E(- 1 ) m_I " l£,m "'" 1 ^W £,i /W- 

^ ''1=0 

The prototypes of spline dQIs are the various univariate and multivariate exten- 
sions of Schoenberg-Marsden operators [52] [53]. 

jeJ 

where tj is an interior point of Sj = supp(B a ). 
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The prototypes of spline iQIs are the various univariate and multivariate extensions 
of operators [21] [63] 

Tf = Y J (LM )B 0l 

jeJ 

where Mj is a B-spline (which can be different from Bj ) normalized by J s Mj = 1 . 

As emphasized by dc Boor ([10], chapter XII), a spline QI defined on non uniform 
partitions has to be uniformly bounded independently of the partition (abbr. UB) 
in order to be interesting for applications. Therefore, with some coworkers, we 
have defined various families of QIs satisfying this property and having an infinite 
norm as small as possible. In general it is difficult to minimize the true norm of the 
operator, however, it is often possible to minimize an upper bound of this norm: 
this gives rise to what we have called near-best (abbr. NB) QIs (see [l],[2]-[4],[40]). 
Numerical applications are still not very much developed. However, QIs can be use- 
ful in approximation and estimation [21] [22] [85] , in numerical quadrature [23] [73] [75] .| 
and for the numerical solution of integral or partial differential equations. 



2 Univariate polynomial QIs 

2.1 Basic operators 

1) The Bcrnstein-Stancu QI [83] is defined for x G [0, 1] by 

n 

s^f(x) = Y J fC-)t ) M 

i=0 

( \ (x) k (1 — x) n ~ k 
where the Bernstein-Stancu basis is defined by oj^'(x, a) = C k a n — - — . 

Here (x) k = x(x + a) . . . (x + (k — l)a),for a G K. For a = 0, we recover the 
classical Bernstein basis. 

2) The Bernstein- Phillips (or q-Bernstein) QI ([56]- [59]) is defined for x G [0, 1] 

by 

^/(*) = E/(|f) tf°(*,ff) 

where the Bernstein-Phillips or q-Bcrnstein basis is defined for q ^ 1 by 
b£\x,q) = - xY q -K Here = W - nUiW, ^ - ^g^jT and 

(x)g = ris=o(l — 1 Sx )- F° r ? = 1, we recover the classical Bernstein basis. 

Using the notation e s (x) = x s for monomials, it is easy to prove that all the above 
QIs B n are exact on Pi, i.e. B n e s = e s for s = 0,1. Moreover they are degree 



4 



preserving since B n e s (x) = e s (x) + r s _i (x, n) where r s -\{x,n) is some polynomial 
of degree at most s — 1 depending on n (and eventually on the parameters a, N 
or q). 
QIs— 

2.2 Left and right BQIs 

All operators B n defined above are isomorphisms of P„. Moreover B n and A n = 
B~ x can be expressed as linear differential operators with polynomial coefficients 

n n 

B n = Y,P { k n) D\ A n = J2^D k , 

k=0 k=0 

where D = -jj^ and the polynomials £ Pfe and 4 n) e P fc are defined by simple 
recursions (see e.g.[64]-[67] for partial results in this sense). 
For < r < n, we introduce the partial inverses: 

fc=0 

and we consider the two families of right and left BQIs: 

(RBQI) The right BQIs B l n ] =B n o arc defined for C r -functions / by 

fiM/ = B n {A[:\f) - B n (£a k n) D k f). 

k=0 

(LBQI) The left BQIs B ( n ] = A^ ] oB n are defined on any (e.g. continuous) function 

r 

Btff = A^(B n f) = Y,a[ n) D k {B n f). 

\r] (r) 

By construction, for < r < n, the BQIs Bn and B„ are exact on the space 
P r . Moreover, in many cases, the LBQIs have a uniformly bounded infinite norm, 
independent on n for each < k < n fixed (see e.g. [70] [86] for some results of 
this type). From this property are deduced some convergence results (see [30] [67]). 



2.3 Kageyama QIs 

Kageyama [44] [45] considers Stancu operators for a e [— -,0] 
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where C n is the Lagrange interpolation operator on the uniform partition of [0, 1] 
(this result is due to Muhlbach). Then he truncates at order s the Maclaurin series 
of Sk a) f w.r.t. a and he takes the value of this polynomial at a = —-: 



4-) / = £ * H£ ^ [5(a)/ 

j=o J 



a=0 



K,n^ = Sn' 1 = B n and /cl°°' 1 = Sn = £„. He also gives expansions of K n f in 
terms of derivatives of B n f and in powers of K He proves that, for all s fixed, 

ll^n^lloo is uniformly bounded and give Voronovskaja type results, e.g. 

^n^(/CW/-/) = -E^^ fe / 

fe=0 

where the polynomials Ts+i^ can be computed by recursion. He also compares 
the expansions of £„, the BQIs and ICn^ interms of derivatives of B n f with 
polynomial coefficients. Numerical experiments done by the author suggest that 
these operators are in general better approximants than BQIs of section 3.1. 



2.4 Univariate Durrmeyer and Goodman-Sharma QIs 

A straightforward generalization of the Durrmeyer operator M n consists in intro- 
ducing a Jacobi weight on [0, 1] in the associated scalar product 

(f,g)= f w a ,f,(t)f(t)g(t)dt, w a , (t)=t a (l-t) 13 , for a,0>-l 
Jo 

The extended Durrmeyer- Jacobi operator ([6], [61]) is then defined by 

1 



The limit case (a, (3) = (—1,-1), corresponding to the weight w(x) = 

x(l — x) 

gives a QI with very attractive properties. It has been introduced by Goodman 
and Sharma [38] [39] for polynomial (and a variant for spline) QIs. It can be written 
as follows, with Lf(x) = (1 - x)f(0) + xf(l): 

n-l 

G n f = Lf+(n-l)J2(f-L.f,bti 2) )t ) - 
»=i 

This operator is exact on Pi and its behaviour is quite similar to that of the 
classical Bernstein operator. For example, one has for / G C 2 (T) 

limn(/(i) - G n f{x)) = x(l - x)f'(x). 
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It also preserves the positivity, the monotonicity and the convexity of /. As discrete 
Bernstein operators, the above operators G„ have associated QIs in the sense of 
section 3.1 [70] . 

2.5 Extrapolation 

All operators B n described in this section have asymptotic expansions of type 

for some index r, the ip^ (/, x) being linear differential operators depending on 
n and k. Therefore they are good candidates for extrapolation methods (see e.g. 
[14] and [82]). Numerical experiments done by the author show that Richardson 
extrapolation is efficient while the use of variants of epsilon or A 2 algorithms often 
introduce spurious poles in the interval of definition. 



3 Polynomial QIs on a simplex 

3.1 Bernstein operator and associated QIs 

The simplex S of dimension d — 1 is defined in barycentric coordinates as 
S = {x = (xi,x 2 , ■ ■ ■ ,x d ) : |x| = 1} with |x| = £\ =1 \xi\. 

The associated simplex of indices, monomials and partial derivatives are defined 

by E„ = {i = (ii,i 2 , ■ ■ ■ , id) ■ |i| = «■}, X n = {± : i G £„} C S, 

i! = ^ 1 !^ 2 ! . . .,i d \, x ! = x^x% . ..a#, D 1 = D?D% ...£)*- with D s = 

The Bernstein basis of P„ (space of polynomials of total degree at most n) and 

the Bernstein operator are defined respectively by: 

6< n) (x) = ^x ; for i e S„, B n /(x) = ]T /(>^(x) 

As Y^ies n & i" } ( x ) = 1 and E ie s„ ^ & i ( " } ( x ) = x £ = = a; s , for 1 < s < d, where 
e s = (0, 0, . . . , 1, . . . , 0), then B n is exact on Pi. 

Let Z™, i € S n } be the Lagrange basis of P„ associated with the data points A„. 
Then = 5y implies = &["\ hence £?„ is an isomorphism of P„. 

For / e C 2 (5"), we have the Voronovskaja type result ([47] [79] [85]), 

limn[S„/-/] = ^/ 
where 2?/ is the differential operator 

'Df(x) = ^2xiXj(di - dj) 2 . 

i<j 
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B n and its inverse A n = B n 1 in P„ can be expressed as linear differential operators 

ies„ ies„ 

whose coefficients can be computed by recursion. For < k < n, define partial 
inverses 

A<P = £ a^D\ 

As in section 2.2 for univariate QIs, we can consider the two families of operators: 
left Bernstein quasi- interpolants (LBQIs) = o B n , and right Bernstein 
quasi-interpolant (RBQIs) B [ „ ] = B n o A^ ] , where B ( n ] = B l ° ] = B n = and 
B^ = i?lr' — C n — Lagrange interpolation on X n . 

(2) 

We have proved [65] that \\Bn ||oc < 2e? + 1 for all n > 2, and we conjecture that 
for all k > 0, there exists a constant Ck{d) such that for all n > k, 

Halloo < C fc (d). 

We also conjecture the Voronovskaja-type results 

limn r+1 (B< 2r) / - /) = A 2r f, limn r+1 (B^f - f) = A* 2r+1 f, 

where A 2r and ^2r+i are linear differential operators, and the asymptotic expan- 
sions 

Cr+l . C r+ 2 



B^fandB^f~f+^± + 



r+l „r+2 



3.2 Durrmeyer-Jacobi QIs on a simplex 

One can introduce a Jacobi weight on the simplex in the scalar product of L^S): 
w a (x) — x a , (f,g) = J s w a (x)f(x)g(x)dx, and define the Durrmeyer-Jacobi 
quasi-interpolants (DJQIs) 

M f = V {f ^ ]) b (n) 
*G£„ ( e oA ) 

Its eigenvectors are the Jacobi polynomials on the simplex. There holds a 
Voronovskaja type result [13] [79] 

limn(M n f(x) - f(x)) =V a f{x) 

where the differential operator V a is defined by 

V a = x~ a Y^(di - d^Xixfidi -dj). 

i<j 



8 



As M n is an isomorphism of P n , one can expand M n = X)fe=o Sies k @\ D 1 and 

L n = M~ x = X)fc=o SiGS„ oif^D 1 . As in the univariate case [71], the polynomials 

/3j (n) and a\ n ^ are probably linear combinations of Jacobi polynomials on S ([26]) 

Setting L ( n ] = Efe=oE ie s n a\ n) D\ one can define the left DJQIs M^ r) = L { n ] o 

M n , and the right DJQIs M [ n ] = M n o L^, with M^ 0) = M„ and M^ rl) = P n = 
orthogonal projector on P„ in L 2 (S). They have the same properties as univariate 
QIs, and it would be interesting to have detailed proofs, those of [65] [66] being 
only sketched. However, the author thinks that the following operators are still 
more attractive. 



3.3 Jetter-Stockler operators on a triangle 

For the sake of simplicity, we describe them over a triangle (with barycentric 
coordinates {Ai,A2,A3}) in the case of the Legendre weight (w — 1, see [42] for 
the general study on a simplex with Jacobi weight). Using the following notations: 

Aj = dj - di, i < j, D = {D12, 1)13,1)23}, A = {AiA 2 , A1A3, A 2 A 3 }, 

k = (fcia.fcia.M G N 3 , V k = D^D^D^f, 
A k = (A 1 A2) fel2 (A 1 A3) fcl3 (A 2 A3) fc23 , 
the authors define the following basic differential operators: 

1 ™ 

u k = -(-i)i k ^ k A k p k , ^ = -£t/ k , y„ = 2(^)" 1 w/ 

' |k|=£ £=0 

Let M n be the Durrmeyer operator, then they prove that Uu commute with M n 
for all pairs (k, n) and that y n is the inverse of M n in the space of polynomials 
P„. Now, for < r < n fixed, they define partial inverses and left Jetter-Stockler 
quasi-interpolants (LJSQIs) 

r 

y { n ] = ^T(^)- 1 ^, AfW = y£) Mn . 

1=0 

One can also define right JSQIs = M n yn \ Both operators and M^} arc 
exact on P r . Aloreover, for r fixed, the left JSQIs have uniformly bounded infinite 
norms w.r.t. n. Finally, the authors prove Voronovskaja-type results: 

Urn C r n (f-M£-Vf)=U r f 
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3.4 Extrapolation 

All operators B n described in this section have asymptotic expansions of type 

B n f{x) » f{x) + J2 ^H^ 1 
k>r 

(see e.g. [47] and [85]). In particular, the latter reports interesting numerical results 
on Richardson extrapolation of classical Bernstein operators on the triangle. It 
would be interesting to compare these results with those which could be obtained 
by extrapolating the above QIs. 
PARTITIONS 



4 Univariate spline QIs on uniform partitions 

4.1 Univariate differential and discrete QIs 

For the construction of QIs with optimal approximation order, we refer to [15] 
and [16], where general solutions are given, thus completing the initial work by 
Schoenberg in [80] . 

4.2 Near-best spline dQIs 

Consider the family of spline dQIs of order 2m depending on nil arbitrary 
parameters a — (a , ai, . . . , a n ), n > m: 

Q a f = J2 A fd) M 2m(x-i) 

with coefficient functionals 

n 

A/(i) = a f(i) + ]T aj (f(i + j) + f(i - j)) . 

3 = 1 

Setting v(a) = |a | + \ a j\i tnen we nave ||Qa||oo < v{a). By imposing that 

Q a be exact on P r , with < r < 2m — 1, we obtain a set of linear constraints: 
a eV r C M™ +1 . We say that Q* = Q a * is a near best dQI if 

v{a*) = min{^(a);a G V r }. 

There is existence, but in general not unicity, of solutions. 

Example: cubic splines (see [40]). There is a unique optimal solution for n > 2: 

a* = 1 + ^ - < = -g^a ' a*=0 for l<j<n-l 

Moreover, for all n > 4, || Q* 1 1 oo < 1 + 3^2- Here are the first values of ||<5*||oo an d 
v(a*)\ n = 1 : 1.222 & 1.666; n = 2 : 1.139 & 1.166; n = 3 : 1.074 & 1.074. 
PARTITIONS— 
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4.3 Near-best spline iQIs 

A similar study can be done for integral spline QIs. We refer to [2] [40] and we only 
give an example given in these papers. Setting a — (ao,ai, . . . , a n ), n > m and 
Mi(x) = M2m(x—i), we consider Q a f — Xiiez A/(i)Mj with coefficient functional 

n 

A/(i) = a (f, M 4 ) + *o «/> Mi-i) + (/, Mi+j)) . 

As in section 4.2, we have ||<3a||oo < f( a ) and we say that Q* = Q a . is a near best 
iQJif v(a*) = min{j/(a);a G V r }. There is existence, but in general not unicity, of 
solutions. 

Example: cubic splines (sec [40]). There is a unique optimal solution for n > 2: 

2 1 

a ° = 1 + 3^? ' a * n = ~3^2 ' a J = 1 - J ~ U ~ 1 

Moreover, for all n > 4, ||<5*||oo < 1 + 3^7 • Here are the first values of ||<3*||oo an d 
v{a*); n = 1 : 1.5278 & 2.333; n = 2 : 1.2778 & 1.333; n = 3 : 1.1481 & 1.1482. 



5 Bivariate spline dQIs on uniform partitions 

5.1 A general construction of dQIs 

Let ip be any kind of bivariate B-spline on one of the two classical three- or four- 
directional meshes of the plane (e.g. box-splines, see [7], [12], [19]). Let X = supp((p) 
and E* = SflZ 2 . Let a be the hexagonal (or lozenge=rhombus) sequence formed 
by the values {(p(i),i G £*}. The associated central difference operator V is an 
isomorphism of P(v?)> the maximal subspace of "complete " polynomials in the 
space of splines S(<p) generated by the integer translates of the B-splinc <p (see 
[12], [69], [71], [72]). Computing the expansion of a in some basis of the space of 
hexagonal (or lozenge) sequences amounts to expand T> in some basis of central 
difference operators. Then, computing the formal inverse V^ 1 allows to define the 
dQI 

Qf= ^^/(fcM-*) 

fcez 2 

which is exact on F((p). Let us now give two examples which are detailed in [40]. 

5.2 Near-best spline dQIs on a three direction mesh 

Example: let ip be the C 2 quartic box-spline. Let H s be the regular hexagon with 
edges of length s > 1, centered at the origin (here S = iJ 2 ) and let H* = H s f\l?. 
The near-best dQIs have coefficient functional with supports consisting of the 
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center and the 6 vertices of H*, s > 1. The coefficients of values of / at those points 
are respectively 1 + ^ and — , therefore the infinite norm of the optimal dQIs 
Q* is bounded above by v* = 1 + \. Here are the first values of ||<5*||oo and v*; 
n=l: 1.34028 & 2; n = 2 : 1.22917 & 1.25; n = 3 : 1.10185 & 1.111. 

5.3 Near-best spline dQIs on a four direction mesh 

Example: let ip be the C 1 quadratic box-spline. Let A s be the lozenge (rhombus) 
with edges of length s > 1, centered at the origin, and let A* = A s fl Z 2 . The 
near-best dQIs have coefficient functionals with supports consisting of the center 
and the 4 vertices of A*, s > 1. The coefficients of values of / at those points are 
respectively 1 + ^ and — therefore the infinite norm of the optimal dQIs 
Q* is bounded above by v* = 1 + \. Here are the first values of ||<5*||oo and v*\ 
n = 1 : 1.5 & 2; n = 2 : 1.25 & 1.25; n = 3 : 1.111 & 1.111. 



6 Univariate spline QIs on non uniform partitions 
6.1 Uniformly bounded dQIs 

Let us only give an example: we start from a family of DQIs of degree m which 
are exact on P2. 

jeJ 

We recall the expansion [52] [53] 

On the other hand, ^D 2 f(9j) can be replaced on the space P2 by the second order 
divided difference [6j-i,9j,6j+i]f, therefore the dQI defined by 

Qlf = 5>f (f)B jt lif\f) = f(6j) - AfVi-uOjA+Af* 

is also exact on P 2 . Moreover, one can write 

vf\f) = difi-i + hfi + Cifi+i 

with a, = -Af Va^_i(A^_! + Aft), c, = -4 (2) /A^(A^-i + A0 4 ), and 
bi = 1 + ^4- 2 '/A#i_i A#j, So, according to the introduction 

A (2) 

\\Q* 2 \\ o<rr^x(\a i \ + \bi\ + \c i \) <l + 2rnax ' . 
The following theorem [4] extends a result given for quadratic splines in [4] [73] [75] . 
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Theorem 1. For any degree m, the dQIs Q\ are UB. More specifically, for all 
partitions of I: 

Moo<[^(m + 4)] 



6.2 Uniformly bounded iQIs 

General types of integral QIs are studied in [21] [63] [68]. Here, we have chosen to 
study a family of QIs that we call Goodman-Sharma type iQIs, as they first appear 
in [38] . They seem simpler and more interesting than those we have studied in [68] . 
The simpler GS-type IQI can be written as follows 

n+m — 2 

G 1 f = f(t )B + £i(/)£i + /(t n )B„ +m -i, 

i=l 

where the integral coefficient functionals are defined by 

- / Mi_i(t)/(t)dt, 
Jo 

Mj_i(t) being the B-splinc of degree m — 2 with support E»_i = [U- m +i,ti], 
normalized by fi^ — /Uj(eo) = Jq 1 M»_i(t) = 1. It is easy to verify that Gi is exact 
on Pi and that ||Gi||oo = 1- We shall study the family of GS-type iQIs defined by 

n+m— 2 

Gif = f(t Q )B + + b iM.f) + <Hfri+i(f)]Bi + /(t„)-B„ +m _i, 

i=l 

which are exact on P 2 . The three constraints G 2 efe = e^, fc = 0, 1, 2, lead to the 
following system of equations, for 1 < i < n + m — 2: 

(2) (2) (2) (2) 

cii + bi + Ci = 1, Vi-idi + #A + = 0i, p-l^di + jj\ bi + Mi+iCi = ^ • 

This is a consequence of the following facts 

Ai(ei) = / tMi^i(t)dt = —y^U- m+s = 9i, 
Jo rnj^ 

= W(e 2 ) = ^ t'M^dt = m(TO 2 +1) I 2 (T t ) 
2 x ^ 

yr / ^ ti—m+rti—m+s 



m(m +1) ^ 

v 7 Kr 



<s<m 
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Theorem 2. For any degree m, the iQIs G2 are UB. More specifically, for all 
partitions of I: 

||G 2 || 0O <5 
The detailed proof will be given in [78] . 

6.3 Near-best dQIs 

Let us consider the family of dQIs of degree m defined, for the sake of simplicity, 
on I = M. endowed with an arbitrary non-uniform increasing sequence of knots 
T={U;ieZ}, by 

Qf = Q P , q f = J2^ B i- 

Their coefficient functional depend on 2p + 1 parameters, with p > m: 

/*(/) = E Ai(*)/(0 i+a ), 

s — — p 

and they are exact on the space P g , where q < min(m, 2p). The latter condition is 
equivalent to Qe r = e r for all monomials of degrees < r < q. It implies that for 
all indices i, the parameters Aj(s) satisfy the system of q + 1 linear equations: 

£ Ai(s)0[ +s =0l r) , 0<r<q. 

s — —p 

The matrix V t G M (<z+1) x (2p+1) of this system, with coefficients ^(r,s) = 6 r i+a , 
is a Vandermonde matrix of maximal rank q + 1, therefore there are 2p — q free 
parameters. Denoting 6, G the vector in the right hand side, with components 
bi(r) — 6^ r \ < r < q, we consider the sequence of minimization problems, for 
i G Z: 

min||Aj||i, V i \ i = b l . 

We have seen in the introduction that V\{Q) = maxjgzmin ||Aj||i is an upper 
bound of ||Qg||oo which is easier to evaluate than the true norm of the dQI. 

Theorem 3. The above minimization problems have always solutions, which, in 
general, are non unique. 

The objective function being convex and the domains being afhne subspaces, these 
classical optimization problems have always solutions, in general non unique. 
Example of optimal dQIs are given in [1][4][40]. 
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7 Bivariate quadratic spline dQIs on non uniform 
criss-cross triangulations 

At the author's knowledge, the only bivariate box-splines which have been ex- 
tended to non uniform partitions of the plane are C^-quadratic box-splines on 
criss-cross triangulations [20] [62]. Recently, we have constructed a set of B-splincs 
generating the space of quadratic splines on a rectangular domain and we have 
defined a discrete quasi-interpolant which is exact on P2 and uniformly bounded 
independently of the partition [74]- [76]. 
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